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Preface 


This  report  serves  two  purposes.  First,  it  serves  as  the  Users’  Manual  for  the 
SHARC/SAMM  Atmosphere  Generator,  a  stand-alone  program  that  generates 
atmospheric  profiles  for  Air  Force  ER  radiation  codes  that  account  for  systematic 
variability  of  the  atmosphere.  For  this  reason  it  includes  the  results  of  work  done  by 
scientists  at  Spectral  Sciences,  Inc.,  who  incorporated  the  SNOE  database  for  NO  density, 
URAP  H2O  database,  revised  O  atom  profiles ,  and  revised  H  atom  profiles,  in  addition  to 
documenting  the  work  of  Air  Force  Research  Laboratory  scientists,  who  incorporated  the 
NRLMSISE-00  database,  incorporated  choices  for  Geodetic  or  Geocentric  Latitude  input, 
and  tested  and  verified  the  code.  Second,  this  report  also  serves  as  the  final  scientific 
report  for  in-house  work  unit  3352BDA1. 

The  authors  thank  R.  Kennett  (Spectral  Sciences  Inc.)  for  numerous 
improvements  to  the  code,  G.P.  Anderson  (AFRL)  for  providing  the  NRL  Trace 
Climatology  database,  R.  Rodrigo  (Instituto  de  Astrofisica  de  Andalucia)  for  providing 
theoretical  profiles,  and  R.  Sharma  for  valuable  discussions.  The  work  done  by  Spectral 
Sciences,  Inc.  was  funded  under  Air  Force  contract  F19628-93-C-0050. 


1.  INTRODUCTION 


The  Air  Force  Strategic  High-Altitude  Radiance  Code  (SHARC)  and  the  SAMM  (SHARC 
And  MODTRAN  Merged)  code  have  been  developed  to  address  the  strategic  requirements  for 
modeling  infrared  (IR)  background  radiation  and  structure  in  the  upper  atmosphere.  One  of  their 
most  critical  applications  is  in  modeling  radiance  variations,  which  can  occur  over  a  wide  range 
of  spatial  and  temporal  scales.  Short-term  and  small-scale  variations  associated  with  random 
processes  can  be  characterized  and  predicted  statistically.  On  the  other  hand,  systematic 
variations,  which  can  be  predicted  deterministically,  may  be  quite  large  and  thus  play  an 
important  role  in  setting  the  overall  background  radiance  level  for  a  given  bandpass.  In 
particular,  at  the  solar  terminator,  large  radiance  variations  can  occur  over  a  small  (several- 
degree)  range  of  solar  zenith  angle  due  to  photochemical  processes  in  the  atmosphere. 

The  SHARC  Atmosphere  Generator  (SAG)  has  been  designed  to  allow  the  major  known, 
systematic  variabilities  in  the  atmosphere,  including  terminator  and  other  diurnal  effects,  to  be 
practically  incorporated  in  strategic  IR  radiance  calculations.  SAG  is  presently  implemented  as  a 
FORTRAN  subroutine,  which  may  be  run  via  the  supplied  or  user  furnished  driver  program. 
This  allows  SAG  to  be  run  interactively  or  in  a  batch  processing  mode  by  looping  over 
atmospheric  dependencies  in  the  driver  program.  SAG  may  be  used  to  generate  an  atmosphere 
file  for  use  with  MODTRAN  and  a  file  of  species  and  kinetic  temperature  profiles  compatible 
with  SHARC  and  SAMM.  The  profiles  are  customized  for  the  geophysical  and  geographic 
information  input  by  the  user. 

Using  information  on  the  day  of  the  year,  local  time,  solar  activity  indices,  etc.,  SAG 
reasonably  models  the  systematic  variabilities  in  CO2,  O3,  OH,  NO,  H2O,  and  O  atom  densities. 
For  other  species,  diumally-averaged  profiles  are  taken  from  recent  databases.  To  facilitate  use 
without  detailed  inputs,  defaults  are  provided  so  that  simple  designators,  such  as  day/night, 
season,  and  latitude  region  (low,  mid,  or  high),  can  be  specified  as  desired. 

SAG  draws  on  several  existing  empirical  atmosphere  models.  Either  MSISE-90(1)  or 
NRLMSISE-00(2)  may  be  used  for  the  temperature  and  major  species  profiles.  They  provide 
profiles  for  species  including  N2,  O2,  O,  and  H  as  a  function  of  altitude,  latitude,  longitude, 
universal  time  (UT),  local  solar  time  (LST),  daily  average  Ap  index,  and  the  FI  0.7  (previous 
day)  and  F10.7A  (81  day  centered  average)  solar  flux  indices.  The  second  atmosphere  model  is 
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the  NRL  climatology  database(3)  for  altitudes  up  to  120  km.  This  database  is  used  for  the 
SHARC  and  SAMM  species  CH4  and  CO,  and  for  lower  portions  of  the  03  and  O  profiles,  as 
well  as  for  the  additional  species  N2O,  NO2,  and  HNO3  used  in  SAMM.  The  NRL  database^  5 
contains  mean  monthly  mixing  ratios  at  1  to  5  km  altitude  increments  and  10°  latitude 
increments.  SAG  interpolates  between  these  values  and  converts  to  number  densities  using  the 
MSISE-90  or  NRLMSISE-00  total  densities. 

New  water  vapor  climatology  has  been  introduced  into  SAG  2.0.  The  new  climatology  was 
developed  by  the  UARS  Reference  Atmosphere  Project  (URAP),  whose  aim  is  to  provide  a 
comprehensive  reference  description  of  the  stratosphere  based  on  the  data  recorded  by 
instruments  on  the  NASA  Upper  Atmosphere  Research  Satellite  (UARS).(4)  The  URAP  water 
vapor  climatology  was  constructed  from  the  HALOE  (HALogen  Occulation  Experiment/  \ 
MLS  (Microwave  Limb  Sounder  Experiment)/65  and  SAGE  II  (Stratospheric  Aerosol  and  Gas 
Experiment/75  data  and  presented  as  monthly  zonal  means,  where  the  zones  are  designated  by 
latitude  and  pressure/85  The  new  climatology  is  a  substantial  improvement  on  CIRA-1996(95, 
which  is  based  on  pre-URAP  data. 

The  Student  Nitric  Oxide  Explorer  (SNOE)  database  has  been  introduced  to  provide  nitric 
oxide  profiles  between  97  km  and  150  km/l0)  The  database  consists  of  measurements  of  nitric 
oxide  density  in  the  thermosphere  for  the  period  March  11,  1998  to  September  30,  2000.  The 
data  covers  the  latitude  range  of  80°  S  to  80°  N  at  5°  intervals  and  the  longitude  range  180°  W  to 
180°  E  at  24°  intervals. 

The  remaining  species  profiles,  including  those  for  CO2,  OH,  SO2,  and  NH3  are  derived  from 
a  combination  of  standard  concentrations015  used  in  MODTRAN,  photochemical  or  empirical 
models  based  on  inputs  or  outputs  (CO2),  or  some  combination  of  these  (OH),  as  described  in 
this  report. 

The  following  sections  provide  a  detailed  description  of  the  SHARC/SAMM  Atmosphere 
Generator,  including  operating  instructions,  computation  methods,  and  examples.  Appendix  A 
summarizes  the  changes  between  SAG  0.92  and  SAG  2.0,  Appendix  B  summarizes  the  program 
structure,  and  Appendix  C  presents  sample  inputs  and  outputs. 


-2- 


2.  USAGE 


2.1  Entering  Parameters 

Parameters  required  for  the  calculations  are  passed  to  SAG  via  a  driver  program.  In  the 
interactive  mode,  parameters  may  be  entered  by  the  user  during  an  interactive  session  with  SAG. 
A  sample  session  is  illustrated  in  Appendix  C.  When  the  SAG  driver  is  run,  it  first  retrieves  the 
parameters  from  the  previous  session,  which  are  stored  in  the  input  file  “drivstor.dat”.  When  the 
main  menu  appears  (see  C.l  Sample  Interactive  Session),  the  user  responds  by  typing  an  integer. 
Two  of  the  responses  (“3”  and  “4”)  lead  to  sub-menus.  Typing  “0”  or  “0,0”  from  a  sub-menu 
returns  the  user  to  the  main  menu.  New  to  SAG  2.0,  in  the  Main  Menu,  is  the  option  “6  Run  QA 
Tests”.  This  option  will  cause  SAG  to  generate  150  different  atmospheres  and  then  run  a 
comparison  between  these  results  and  150  previously  calculated  atmospheres  in  the  file 
QADAT/ATM_QA.DAT.  The  QA  Test  option  may  be  used  to  test  the  integrity  of  a  SAG 
installation.  When  the  QA  test  is  performed  the  “drivstor.dat”  file  is  overwritten. 

Table  1 :  Parameter  Menu 

PARAMETER  MENU 
Current  parameters  are 

1  Daily  Ap  =  18.0 

2  F10.7  =  120.0 

3  F10.7A  (81  day  centered  average  of  F10.7)  =  120.0 

4  Latitude  Type  =  1  Latitude  (- =S)  =  16.6 
Latitude  Type:  Geocentric(O),  Geodetic(l) 

5  Longitude  =  25.0 

6  UT  (sec)  =  43200. 

7  Solar  Time  (hrs)  =  13.67 

8  Month  =  5 

9  Day  =  12 

10  Year  =  2002 

11  Day  Number  =  132 

12  Solar  Zenith  =  24.76 

13  H,  O,  URAP,  SNOE,  NRLMSISE-00  =  T  T  T  T  F 

Enter  item  number  [n]  or  [0]  to  return  to  main  menu: _ 


The  sub-menus  are  used  to  display  the  current  parameters  and  to  enter  new  ones.  The 
Parameter  Menu  lists  the  current  values  and  permits  the  user  to  directly  enter  new  values.  An 
example  Parameter  Menu  is  shown  in  Table  1.  The  menu  contains  information  on  the  day  of  the 
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year,  universal  time  or  local  solar  time,  latitude,  longitude,  FI 0.7  daily  solar  flux  for  previous 
day  and  F10.7A  81  day  centered  average  solar  flux,  and  daily  average  Ap.  The  SAG  version  2.0 
parameter  menu  also  contains  new  options  for  using  the  NRLMSISE-00  empirical  model  in 
place  of  MSISE-90,  including  the  URAP  or  SNOE  databases,  and  adding  corrections  to  the  O 
atom  and  H  atom  profiles.  The  latitude  may  be  entered  in  either  geocentric  or  geodetic 
coordinates.  SAG  uses  geodetic  coordinates  in  its  calculations. 


Table  2:  Default  Parameter  Values 


PARAMETER 

UMTS 

DESIGNATOR 

VALUE 

Day 

12a 

Local  Solar  Time 

Hours 

Night 

0a 

Winter 

66 

Spring 

157 

Day  of  Year 

Summer 

248 

Fall 

339 

Equatorial 

5 

Latitude 

Degrees  N 

Mid 

35 

High 

65 

Low 

70 

Solar  Flux 

F10.7,  F10.7A 

Medium 

120 

High 

200 

Low 

3 

Geomagnetic  Activity 

Ap 

Medium 

15 

High 

75 

aLongitude  also  set  to  O' 

D 
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The  Default  Menu  is  used  for  entering  default  parameters,  which  are  convenient  for  generic 
calculations,  and  for  specifying  a  custom  list  of  altitudes  (see  Subsection  2.3).  The  default 
values  are  given  in  Table  2.  For  example,  the  geocentric  latitude  is  set  to  35°  when  the  latitude 
default  designator  “m”  (mid-latitude)  is  entered.  The  letter  “m”  will  then  appear  at  the  end  of  the 
“latitude”  line  in  the  menu  listing.  The  absence  of  a  designator  letter  indicates  that  a  directly- 
entered  parameter  is  in  use. 

One  can  re-enter  any  parameter  or  either  sub-menu  any  number  of  times.  By  using  both  sub¬ 
menus,  any  combination  of  default  parameters  and  directly  entered  parameters  can  be  chosen. 
To  save  a  particular  set  of  parameter  values  for  future  re-use,  one  can  save  a  copy  of 
“drivstor.dat”  under  a  different  file  name.  This  file  can  then  be  recopied  to  “drivstor.dat”  when 
the  SAG  rerun  is  needed. 

2.2  Parameter  Relationships  and  Permitted  Values 

As  there  are  more  parameters  in  SAG  than  there  are  independent  variables,  some  parameters 
are  constrained.  For  example,  Universal  Time  (UT)  and  Local  Solar  Time  (LST)  are  mutually 
dependent  for  a  given  longitude.  When  either  UT  or  LST  are  entered,  the  program  automatically 
updates  the  other  parameter.  When  longitude  is  entered,  LST  is  updated  using  the  existing  value 
of  UT.  The  longitude  range  is  -180°  to  360°  to  allow  the  user  to  follow  either  the  0  to  360°  or 
-180  to  180°  convention,  where  0°  to  180°  is  defined  to  be  West  longitude.  The  parameter 
values  are  also  checked  for  physical  reasonableness. 

Except  when  directly  entered,  the  solar  zenith  angle  (SZA)  is  treated  as  a  dependent  variable, 
being  a  function  of  the  longitude,  latitude,  Julian  day,  and  either  LST  or  UT.  The  program 
automatically  calculates  and  updates  SZA  when  any  of  those  parameters  are  entered.  The  SZA 
computed  by  SAG  is  not  transferred  directly  to  SHARC  or  SAMM  at  present.  When  the  SZA  is 
directly  entered,  the  other  parameters  are  retained  except  for  longitude,  which  is  recomputed. 

2.3  Altitude  and  Species  Lists 

The  list  of  desired  altitudes  and  the  choice  of  SHARC  or  SAMM/MODTRAN  species  are 
specified  in  the  Default  Menu.  Default  lists  exist  for  SAMM,  SHARC,  and  MODTRAN. 
Alternatively,  a  custom  altitude  list  can  be  input  by  the  user  from  file  “drivht.dat”.  Line  1 
contains  an  integer  argument,  0  for  outputting  SHARC  species  only  or,  1  for  outputting  all 
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SAMM/MODTRAN  species  (see  Appendix  B.5).  The  number  of  altitudes  appears  on  line  2.  Up 
to  200  altitudes  from  0  to  300  km  are  allowed.  The  altitude  list  starts  on  line  3;  the  FORTRAN 
general  input  format  is  used,  so  that  either  integer  or  floating  point  numbers  may  be  entered. 
Note  that  the  listing  in  the  output  atmosphere  file  is  rounded  to  two  decimal  places. 


Figure  1 :  A  flowchart  explaining  the  new  SAG  2.0  capabilities  accessible  through  option  13  of 
the  PARAMETER  MENU.  Choosing  F  (false)  for  all  options  returns  the  SAG  0.92  atmosphere. 


2.4  SAG  2.0  Options 

SAG  2.0  has  several  new  capabilities  accessed  through  option  13  of  the  PARAMETER 
MENU.  Figure  1  shows  a  flowchart,  which  describes  many  of  the  new  feature  of  SAG  and  their 
relation  to  each  other.  The  NRLMSISE-00  option  is  used  to  choose  between  the  two  MSISE 
databases,  which  provide  temperature  and  major  species  profiles.  These  databases  along  with 
the  NRL  database(3)  are  used  by  SAG  to  generate  the  minor  species  profiles,  as  described  later  in 
this  manual.  The  SNOE  and  URAP  database  profiles  for  NO  and  H20  may  be  used  over  their 
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defined  altitude  ranges  to  replace  the  SAG  0.92  style  profiles.  Correction  to  the  O  and  H  density 
profile  may  also  be  included,  the  details  of  which  will  be  presented  later  in  the  manual . 


3.  PROFILE  COMPUTATIONS 

This  section  describes  how  the  atmospheric  profiles  are  generated.  The  Temperature,  N2  and 
O2  profiles  are  taken  directly  from  MSISE-90  or  NRLMSISE-00  and  have  been  discussed  in 
detail  in  Hedin's(1,2^  papers.  Brief  discussions  are  given  in  the  following  subsections  on  the 
profiles  for  temperature  and  the  remaining  species. 

Diurnal  variations  in  the  CO2,  NO,  O3,  O,  and  OH  densities  are  handled  in  an  empirical, 
analytical  manner  in  SAG  using  an  interpolation  weighting  factor,  W,  in  combination  with 
limiting  daytime  and  nighttime  profiles: 

X=  WXd  +  (l-W)Xn  (1) 

Here  X  is  the  density  of  species,  Xd  is  the  daytime  (photochemical  equilibrium)  density  and  Xn  is 
the  nighttime  (photochemical  equilibrium)  density  (see  subsequent  sections  for  detailed 
descriptions  for  individual  species).  The  daytime  and  nighttime  densities  are  functions  of  the 
geophysical  and  geographic  parameters  (excluding  time)  and  the  dependencies  of  W  are 
discussed  below. 


For  O3,  O,  and  OH,  the  chemical  time  constants  are  often  comparable  to  or  shorter  than  the 
sunrise  or  sunset  duration.  The  weighting  factor  for  these  species  is  given  by. 


W(0,)  = 


<9-0 


n  __  n 
^max  ^min 

0 


if  Qs  <  0min 

if0min<05<0r 
if  0S  >  0max 


(2) 


where  6S  is  the  solar  zenith  angle  calculated  4.5  minutes  prior  to  the  current  time  in  order  to 
approximate  the  photochemical  lag  time.  0mjn  and  0max  are  the  minimum  and  maximum  solar 
zenith  angles,  which  define  the  range  of  the  day/night  switch.  The  minimum  and  maximum  solar 
zenith  angles  depend  on  the  altitude. 

0mm  =  [90.57 +  0.123(fl/f- 50)]- 2.5  and  0max  =  [90.57  +  0.123(a/f -50)]+2.5  (3) 
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where  alt  is  the  altitude  (km)  and  the  terms  in  the  square  brackets  are  linear  fits  to  the  50%  point 
solar  zenith  angle.  The  factor  of  2.5  provides  a  5  degree  range  and  the  factor  of  0.123  causes  the 
midpoint  of  the  interpolation  region  to  vary  at  a  rate  of  0.123°  per  km.  W  approaches  or  equals  0 
during  the  night  and  1  during  the  day,  and  is  set  to  0  during  continuous  polar  nighttime  and  1 
during  continuous  polar  daytime. 

For  C02  and  NO,  W  is  calculated  with  the  assumptions  that  the  approach  to  the  daytime  or 
nighttime  equilibrium  starts  immediately  at  sunrise  or  sunset  (SZA  =  90°)  and  proceeds  with  a 
characteristic  time  constant,  x.  An  analytical  solution  to  W(7RjSe,  Tset,  Tlst)  exists  in  this  case, 


W  (TRisc,TSel,TLST)  = 


1  +  (wr  - 1  y;,(r*.-r«r )  if  Trlk  <  tlst  <  4 

Wy;‘(r&,-rm)  if  Tlst  >  TSe, 
if  Tlst  <  TRise 


where  Ws  and  Wr  are  the  values  at  sunset  and  sunrise,  Tlst  is  the  local  standard  time,  T Rise  is  the 
sunrise  time,  TSe t  is  the  sunset  time,  7n  is  the  night  time  constant,  and  rd  is  the  day  time  constant 
(7d  =  7n  =  3  for  NO  and  7d  =  7n  =  8  for  C02).  Here, 


0  _  &  ij'sct  1. Rise)'Td 

—  &  , 


1  \}~eden) 

e  _e-(24-(rs,,-rtec))/r„ 


In  calculating  the  sunrise  and  sunset  LST’s,  0°  longitude  is  assumed.  This  can  lead  to  a  small 
error,  typically  less  than  one  minute,  compared  to  the  sunrise  and  sunset  times  at  the  actual 
longitude. 


3.1  SHARC  Species 

3.1.1  C02 

The  C02  profile  is  computed  from  specified  nighttime  and  daytime  ratios  relative  to  the 
argon  number  density,  which  is  obtained  from  MSISE-90/NRLMSISE-00.  Below  90  km,  a  fixed 
C02/Ar  ratio  is  used,  based  on  a  ground  level  C02  mixing  ratio  of  340  ppm.  Above  90  km,  the 
ratio  decreases,  and  separate  daytime  and  nighttime  equilibrium  ratios  are  used.  These  daytime 
and  nighttime  C02/Ar  ratios  were  obtained  by  combining  the  theoretical  calculations  of  Rodrigo 
et  al.(X2)  for  sunrise  (representing  nighttime  equilibrium)  and  sunset  (representing  daytime 
equilibrium)  with  mass  spectrometric  measurements.  MSISE-90  was  used  to  compute  Ar 
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densities  corresponding  to  the  conditions  of  Rodrigo  et  al.(n)\  using  their  CO2  densities,  the 
CCVAr  ratios  were  taken.  These  ratios  were  then  systematically  rescaled  for  overall  consistency 
with  measurements  from  the  daytime  Aladdin  experiment 13)  and  lower-altitude  nighttime 
measurements. (14)  For  CO2  diumally  dependent  interpolating  factors  are  defined  at  120  km  and 
150  km, 


density  ratios. 


X120  =  LOG(0.013W  +  0.18(1  -  W»  and  Xi50=  LOG(0.010W  +  0.17(1  -  W))  (6) 

where  X/,,  is  the  CCVAr  ratio  at  an  altitude  of  ht  km.  The  value  of  W  is  determined  from 
equation  4  with  a  relaxation  time  constant  x  of  8  hrs  for  both  day  and  night.  The  CO2  density  is 
then  given  by, 


[C02]  =  [Ar]X90 

[C02]  =  [  Ar]<?(  ( xl2°(aU-9°) +  X9°(120  -  alt)  V3°) 

[C02]  =  [Ar]e((X150(a"_I20)  +  xl20(150_all))/30) 
[CO2]  =  [Ar]e(  ( X30°(alt_15°) +  X150(300  ~  ait)  )/15o) 


for  alt  <  90  km 

for  90  km  <alt  <  1 20  km 

for  120  km  <alt  <  150  km  (7) 

for  150  km  <alt 
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The  resulting  daytime  and  nighttime  ratios  are  shown  in  Figure  2  along  with  the  theoretical 
calculations  and  the  Aladdin  measurements.  C02  concentration  profiles  generated  by  this  model, 
such  as  the  examples  in  Figure  3,  are  found  to  be  consistent  with  measurements  summarized  in  a 
recent  review/1 5^  and  are  very  similar  to  mid-latitude  profiles  used  previously  for  IR  radiation 
modeling/15, 16)  New  C02  concentrations  have  been  derived  from  ATMOS  data(17)  at  altitudes  of 
up  to  120  km.  The  model  is  found  to  provide  an  excellent  representation  of  the  average  of  these 
somewhat  scattered  measurements.  From  these  comparisons,  the  SAG  C02  profile  is  estimated 
to  be  accurate  to  within  ±30%  or  better  in  the  100-150  km  range. 


Figure  3:  The  SAG  2.0  C02  concentration  profiles  for  autumn,  mid-latitude,  and  moderate- 
activity  conditions. 


3.1.2  NO 

Updates  to  the  NO  profile  based  on  the  SNOE  database  may  be  included  from  the  Parameter 
Menu  (see  Table  1).  When  the  SNOE  database  is  not  used  then  the  NO  profiles  from  SAG 
version  0.92  are  obtained. 
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SAG  0.92  NO  Profiles 

Below  90  km,  the  NO  concentration  is  taken  from  the  AFGL  Constituent  Profiles 
tabulation/1 While  the  actual  concentration  should  have  strong  diurnal  and  latitude  variations, 
little  or  no  NO  radiation  is  expected  to  be  observable  from  this  altitude  range  due  to  competing 
sources.  At  higher  altitudes,  NO  emission  is  quite  strong,  so  the  variations  need  to  be  modeled. 
Above  100  km,  separate  daytime  and  nighttime  profiles  are  calculated  from  the  empirical 
latitude,  temperature,  and  Kp-dependent  models  of  Smith  et  al./18^  as  discussed  below.  The  time 


Figure  4:  A  plot  of  the  mid  latitude  daytime  NO  profile  showing  the  density  interpolation.  The 
altitudes  90  km  and  100  km  are  marked  with  dotted  lines. 

constant  t  is  taken  as  3  hours.  The  densities  between  90  and  100  km  are  estimated  by 
interpolation  (Figure  4).  Therefore,  for  90  km  and  higher  the  NO  density  is  given  by, 

[NO]  =  [N09o]01(100"alt)[NO,oo]  0  -°-iCioo-ait))  for  gQ  ^  ^  <100  ^  (8) 

[NO]  =  W[NOday]  +  (1  -  W) [NOnight]  for  1 00  km  <  alt  (9) 
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where  W  is  from  equation  4,  [NO90]  is  the  NO  number  density  at  90  km  calculated  from  the 
MODTRAN  database,  [NOioo]  is  the  NO  number  density  at  100  km,  [NOday]  is  the  daytime  NO 
profile,  and  [NOnight]  is  the  nighttime  profile.  The  daytime  and  nighttime  profiles  were  chosen  to 
give  consistency  with  the  latitude-dependent  diurnal  variation  in  NO  reported  by  Smith  et  al.  and 
with  low-latitude  satellite  measurements0 9)  during  daylight  hours.  The  nighttime  profile  above 
100  km  is  specified  using  a  fitting  parameter,  the  so-called  equatorial  night/day  ratio 
NOe(n)/NOe(d),  which  is  given  an  altitude-dependent  value  of  between  0  and  0.2.  This  ratio  is  a 
fictitious  quantity  since  the  minimum  latitude  absolute  value  for  the  NO  computation  is  taken  to 


Figure  5:  A  plot  of  the  nighttime  NO  density  as  a  function  of  latitude  for  SAG  0.92  at  1 10  km. 

For  latitude  absolute  values  >  70°  the  latitude  is  “folded  back”  to  give  better 
agreement  with  satellite  data. 


be  10°,  to  ensure  that  finite  values  of  nighttime  NO  are  obtained  at  the  equator.  For  latitude 
absolute  values  >  70°,  the  latitude  is  “folded  back”  (i.e.  effective  latitude  =  140-lat,  see  Figure  5) 
to  give  better  agreement  with  satellite  data  compiled  by  Barth/  ^  Taking  average  nighttime  and 


-  12- 


daytime  NO  densities  as  represented  by  LST  =  0  and  12  hr  respectively,  the  results  at  1 10  km  are 
in  reasonable  agreement  with  the  data  compiled  by  Smith  et  al.(18} 

The  time-dependent  behavior  of  NO  according  to  this  empirical  model  is  shown  in  Figure  6. 
The  low-latitude  simulation  reasonably  reproduces  the  increase  in  mesospheric  NO  at  sunrise 
found  in  the  Atmospheric  Explorer  C  measurements/19'  However,  it  predicts  the  maximum  to 
occur  at  sunset,  whereas  the  observed  maximum  is  somewhat  earlier,  in  the  mid-afternoon.  The 
nighttime  behavior  of  NO  at  low  latitudes  is  not  known  from  measurements^  as.  din.  satellite  data 
are  from  a  solar  fluorescence  technique  and  there  have  been  very  few  rocket  measurements  at 
low  latitudes.  Available  IR  measurements (i8)  show  little  diurnal  variation  at  polar  latitudes. 


SAG  2.0  NO  Profiles 

If  the  SNOE  database  option  is  chosen  then  an  NO  density  profile  based  on  the  SNOE 
database  is  smoothly  integrated  into  the  existing  SAG  profile  for  altitudes  between  97  km  and 
150  km.  The  SNOE  data  were  averaged  to  produce  a  database  that  gives  the  NO  density  as  a 
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function  of  altitude  and  latitude  for  each  month.  The  density  dependence  on  Ap  index  is 
determined  by  correlating  the  SNOE  density  against  an  Ap  index  database.(21)  Variation  in 
density  as  a  function  of  Ap  index  was  fit  via  linear  least  squares  at  specific  altitudes.  These  fits 
were  then  used  to  generate  a  scaling  factor  for  the  average  NO  density,  which  is  a  function  of 
altitude  and  Ap  index  (see  Figure  7).  The  magnitude  of  the  diurnal  variation  is  taken  from  the 
SAG  0.92  profiles  and  used  to  scale  the  SNOE  data.  A  comparison  of  the  NO  profiles  for  SAG 
with  and  without  the  SNOE  database  is  shown  in  Figure  8. 


Figure  7:  The  NO  scaling  factor  for  an  Ap  index  of  15. 


3.1.3  O  Atoms 

An  option  to  scale  and  augment  the  O  atom  profiles  to  improve  agreement  with  experimental 
measurements  has  been  added  to  the  Parameter  Menu  (see  Table  1).  These  updates  may  be  used 
with  either  MSISE  model.  When  MSISE-90  is  used  without  the  updates,  the  O  atom  profiles 
from  SAG  version  0.92  are  obtained. 
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Figure  8:  Nitric  oxide  number  density  differences  between  SAG  0.92  and  SAG  2.0  using 
NRLMSISE-00  and  SNOE  for  noon,  spring,  mid-latitude,  and  moderate-activity 
conditions.  The  figure  shows  the  relative  difference:  ([NO]o.92  -  [NO]2.o)/[NO]2.o, 
where  [NO]  is  the  nitric  oxide  number  density  and  the  subscript  indicates  the  SAG 
version.  The  differences  shown  below  97  km  are  due  to  using  different  MSISE 
models  when  converting  from  mixing  ratios. 


SAG  0.92  O  Atom  Profiles 

The  O  atom  profiles  are  adopted  from  MSISE  above  86  km  altitude.  However,  at  lower 
altitudes  they  are  unrealistic,  as  they  lack  a  diurnal  variation.  Alternative  approaches  were 
therefore  used  for  the  lower  altitudes. 

From  72  to  86  km,  daytime  and  nighttime  profile  shapes,  scaled  to  the  value  at  86  km,  were 
estimated  from  various  measurements  and  calculations  in  the  literature.  The  nighttime  profile  is 
essentially  equivalent  to  that  of  the  U.S.  Standard  Atmosphere.  Below  72  km,  the  daytime  O 
atom  profile  is  computed  from  the  daily  average  value  from  the  NRL  climatology  database'31 
under  the  assumption  that  the  nighttime  value  is  zero;  the  day  length  is  accounted  for.  The  result 
is  rescaled  to  match  the  lower  boundary  of  the  72-86  km  region  daytime  profile.  The  scale 
factor,  which  is  typically  between  0.6  and  1 .6,  is  reported  to  the  user.  The  O  atom  profiles  are 
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switched  from  day  to  night  using  the  linear  intexpolation  scheme  of  the  ozone  terminator  model, 
discussed  in  the  following  subsection.  In  particular: 


[0]  =  W[0]day  +  (l-W)S0night[0]86 
[O]  =  W  x  S0day[0]86  +  (1-  W)S0night[0]86 


[0]day  —  t^NRL 


[0]p72,MSIS  12 

[0]d72  ,NRL  TRise 


alt  =  72  km 
72  km  <  alt  <  86  km 


(10) 


where  the  SO  parameters  are  scaling  factors  for  the  shape  of  the  O  profile. 

The  predicted  O  atom  profiles  should  be  reasonable  at  most  altitudes  during  the  day,  when 
the  concentration  is  nearly  in  steady  state,  but  almost  meaningless  for  altitudes  below  80  km  at 
night,  when  the  concentration  drops  rapidly.  Fortunately,  such  low  O  atom  concentrations  are 
believed  to  have  little  impact  on  atmospheric  IR  radiation. 


SAG  2.0  O  Atom  Profiles 

An  updated  O  atom  profile  may  be  derived  from  either  of  the  MSISE  models.  The  O  atom 
profile  was  empirically  modified  to  be  in  better  agreement  with  experimental  measurements^2' 
28>32),  while  retaining  diurnal  variability.  This  was  done  by  slightly  reducing  the  MSISE  density 
by  scaling  with  a  non-linear  function  of  altitude  near  the  peak  density. 

[O]  =  [0]ws/5£  (5  +  (1  -  6')(abs(alt  -  altmax )/  rf)2 )  for  altmax  -  d  <  alt  <  altmax  +  d  (11) 

where  [0]msise  is  the  MSISE  oxygen  density,  5  =  0.875 -C(l-£>„)  is  a  scaling  factor,  C  =  0.2 
for  MSISE-90  and  C  =  0.3  for  NRLMSISE-00,  altmax  is  the  altitude  at  the  MSISE  maximum 
oxygen  density,  Dn  controls  the  diurnal  variability,  and  d  (=  25  km)  controls  the  scaling  width. 
The  density  is  augmented  by  adding  the  lower  peak  to  better  reproduce  the  O  atom  profile  and 
diurnal  variability, 

[O]  =  [0]lliax  Dne^ ^  2)((a,,~(aU">«"al  W  ))/(FWHM  )  for  ajt  <  altmax  (12) 

where  [0]max  is  the  maximum  in  the  MSISE  density  profile  located  at  an  altitude  altmax,  FWHM 
is  the  full  width  at  half  max  of  the  distribution,  and  a\tpe„k  is  the  altitude  of  the  secondary 
maximum  in  the  oxygen  distribution.  For  Oxygen,  the  FWHM  employed  is  22  km  for  altitudes 
above  the  secondary  maximum  and  34  km  for  lower  altitudes.  Dn  varies  from  1.0  to  0  (day  to 
night).  Comparisons  of  the  original  and  updated  densities  are  shown  in  Figure  9  and  Figure  10. 
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Figure  9:  A  comparison  of  the  original  and  updated  O  atom  density  profiles  for  mid  latitude 
summer. 
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Figure  10:  Oxygen  number  density  differences  between  corrected  and  uncorrected  profiles  using 
MSISE-90  for  spring,  mid-latitude,  and  moderate-activity  conditions.  The  figure 
shows  the  relative  difference:  ([0]corrected  -  [0]uncorrected)/[0]uncorrected, 
where  [O]  is  the  oxygen  number  density. 


3.1.4  03 

The  ozone  profile  is  presently  derived  from  a  combination  of  the  NRL  database(3)  and  a 
photochemical  equilibrium  model.  The  NRL  database(3)  contains  daytime  mixing  ratios  and 
night/day  ratios,  and  is  used  below  72  km.  The  latitude  range  is  restricted  to  ±60°  to  avoid 
difficulties  during  continuous  polar  day  and  night.  The  photochemical  model  is  used  at  80  km 
and  above  for  MSISE-90  and  82  km  and  above  for  NRLMSISE-00.  The  region  between  72  and 
80  km  (or  82  km)  is  connected  by  interpolation. 

SAG  calculates  ozone  density  above  80  km  (or  82  km)  using  the  following  photochemical 
equilibrium  model: 

[03]  =  {ki[0][02]2(300/T)2  +  k2[N2]  [0][02]  (300/T)28+k3[0]2[02]exp(345/T)  +  k4[0][02]} 

/  {k5[0]exp(-2060/T)  +  k6[H]exp(-480/T)},  (13) 

where  for  the  production  (numerator)  terms,  k]  =  6.2x10 34  cm6/(molec2-sec),  k2  =  5.7X10 
cm6/(molec2 •  sec),  k3  -  2.2x1  O'34  cm6/(molec2-sec),  and  U  =  8.0x1  O'21  cm3/(molec-sec);  and  for 
the  loss  terms  (denominator),  k3  =  8.0x10  12  cm3/(molec-sec)  and  kg  =  1.4X10  cm  /(mo  lee  sec). 
The  photochemical  equilibrium  model  accounts  for  ozone  formation  from  O  atoms  and 
destruction  via  H  atom  and  O  atom  reactions  (the  former  yielding  OH)  and  solar  photolysis. 
Standard  rate  constants  from  the  literature  ' :  have  been  adopted.  The  inputs  are  the 

temperature,  N2,  02,  O,  and  H  profiles  and,  if  daytime,  the  photolysis  rate.  As  discussed  in  the 
following  subsection,  uncertainty  in  the  H  atom  profile  leads  to  a  comparable  uncertainty  in  the 
nighttime  ozone.  A  two-body  radiative  association  term  is  included  as  a  formation  mechanism, 
and  has  a  significant  and  increasing  impact  above  100  km.  The  value  of  the  rate  constant, 
8xl0~21  cm3molec‘1s"1,  is  estimated  from  rate  constants  for  formation  and  relaxation  of  the  02-0 
complex(30)  and  from  v3  band  Einstein  coefficients.  It  is  almost  certainly  an  upper  limit. 
However,  for  systems  applications  it  is  safer  to  overestimate  the  ozone  than  to  underestimate  it. 
Furthermore,  any  overestimation  from  this  2-body  channel  would  be  offset  by  the  presence  of 
certain  02*  +  02  reactions  that  have  been  proposed  to  form  ozone(29)  but  which  are  not  presently 
included  in  SAG. 
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Figure  1 1 :  A  comparison  of  the  SAG  model  and  NRL  climatology  ozone  concentrations  for  fall 
equinox,  45°  N,  and  moderate-activity  conditions.  Plots  for  SAG  2.0  include 
corrections  to  H  atom  and  O  atom  densities. 


A  comparison  of  typical  SAG  and  NRL  ozone  concentration  profiles  is  provided  in  Figure 
1 1 .  There  is  quite  good  agreement  between  the  two,  except  for  the  nighttime  profile  for  SAG 
0.92  at  80  km.  This  can  probably  be  ascribed  to  a  slight  underestimation  of  the  H  atom  density 
at  that  altitude  and  by  an  overestimation  of  O  atom  density  in  SAG  0.92.  The  new  O  and  H 
profiles  offered  by  SAG  2.0  (based  on  either  MSISE  model)  agree  much  better  with  the  NRL 
values  at  80  km  but  are  consistently  lower  at  95  km,  which  brings  the  subsequent  predicted 
radiances  into  better  agreement  with  the  CIRRIS  1A  experiment/31^  A  comparison  of  the 
calculated  ozone  density  for  both  MSISE  models  is  shown  in  Figure  12. 
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Figure  12:  Ozone  number  density  differences  for  SAG  0.92  and  SAG  2.0  with  and  without  H 
and  O  corrections  for  spring,  mid-latitude,  and  moderate-activity  conditions.  The 
figure  shows  the  relative  difference:  ([03]COn—  [03]Uncorr)/[03]uncorr,  where  [O3]  is  the 
ozone  number  density  and  the  subscript  indicates  the  corrected  and  uncorrected 
profiles. 


The  transition  between  the  day  and  night  profiles  is  handled  using  a  simple  empirical  model 

/ 1  9  \ 

that  approximates  the  terminator  behavior  in  the  time-dependent  calculations  of  Rodrigo  et  al., 
which  were  performed  for  mid-latitude  equinox  conditions  above  60  km.  The  transition  is  made 
by  linearly  interpolating  over  a  5°  range  of  SZA  computed  for  LST-0.075  hr;  the  0.075  hr  (4.5 
min)  difference  accounts  for  the  approximate  photochemical  lag  time.  The  midpoint  of  the 
interpolation  region  is  taken  to  vary  linearly  with  altitude  at  the  rate  of  0.12°  per  km.  Based  on  a 
comparison  with  Rodrigo  et  al.’s  calculations,  this  simple  representation  of  the  terminator 
location  and  width  should  be  accurate  to  within  around  1°  in  SZA  in  the  70-100  km  altitude 
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region.  Typical  time-dependent  profiles  generated  by  SAG  with  MSISE-90  for  the  dusk 
terminator  are  shown  in  Figure  13. 
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Figure  13:  The  time-dependence  of  ozone  concentrations  at  the  dusk  terminator  for  fall  equinox, 
45°  N,  and  moderate-activity  conditions. 


3.1.5  OH 

The  OH  profile  is  taken  from  a  standard  AFGL  profile  tabulation^1  below  72  km,  a  chemical 
equilibrium  model  above  80  km,  and  interpolated  values  in  between: 


WT  =  (alt  -  72)/8 

(14) 

[OH]  =  [OH72](,0"WT)[OH8o]WT 

for  72  km  <  alt  <  80  km 

(15) 

rom  ki[°3][H]exP(-480/T) 

for  80  km  <alt 

(16) 

[0]exp(-l  17/T)/k2 

where  kj  =  1.4  *  10-10,  k2  =  2.2  x  10” n,  [OH72]  is  the  OH  number  density  at  72  km  and  [OHgo]  is 
the  OH  number  density  at  80  km.  The  chemical  model,  involving  formation  via  the  H  +  O3 
reaction  and  destruction  via  O  atoms,  insures  that  the  H,  O3,  and  O  concentrations  are  self- 
consistent  and  thus  accounts  for  diurnal  variation.  The  standard  profile  neglects  diurnal 
variation.  Little  or  no  sensitivity  of  atmospheric  radiance  to  the  ground  state  OH  density  is 
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expected,  since  it  neither  contributes  to  nor  significantly  absorbs  the  dominant  chemiluminescent 
(Meinel  band)  radiation  from  the  higher  vibrational  states. 


Figure  14:  A  plot  of  the  mid  latitude  daytime  OH  profile  showing  the  density  interpolation.  The 
altitudes  72  km  and  80  km  are  marked  with  dotted  lines. 

3.1.6  H  Atoms 

An  option  to  scale  and  augment  the  H  atom  profiles  based  on  comparisons  with  experimental 
measurements  has  been  added  to  the  Parameter  Menu  (see  Table  1).  When  these  updates  are  not 
used  then  the  H  atom  profiles  from  SAG  version  0.92  are  obtained. 

SAG  0.92  H Atom  Profiles 

The  H  atom  profile  is  taken  from  MSISE-90.  As  with  the  O  atom  profile,  the  H  atom  profile 
cuts  off  below  74  km,  and  its  accuracy  is  suspect  below  80  km  due  to  the  lack  of  diurnal 
variation.  Compared  with  the  nighttime  H  atom  profiles  derived  by  Thomas*32^  for  the  80-93  km 
altitude  range,  the  corresponding  MSISE-90  predictions  appear  to  average  a  factor  of  2  to  4 
lower.  However,  the  uncertainty  in  Thomas’  values  is  reportedly  a  factor  of  2,  and  other 
experiments  ^33,34)  show  better  agreement  with  the  MSISE-90  profile. 

If  the  H  atom  profile  was  indeed  too  low,  the  photochemical  model  for  ozone  and  OH, 
implemented  at  and  above  80  km,  would  predict  too  much  O3  at  night  and  too  little  OH  during 
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the  day.  There  is  some  evidence  for  the  first  effect  at  80  km  (see  Subsection  3.1.4).  However, 
the  daytime  O3  and  nighttime  OH  densities  are  not  affected  by  the  H  atom  density. 

SAG 2.0  H Atom  Profiles 

An  updated  H  atom  profile  may  be  derived  from  either  of  the  MSISE  models.  The  H  atom 
profile  was  empirically  modified  to  be  in  better  agreement  with  experimental  measurements/35' 
41>32)  ^  Gaussian  distribution  of  H  density,  centered  at  the  H  peak  density,  is  added  to  the 

profile. 

[H]  =  [H]nlax  e(-41"(2)((alt-al«mJ/(Z»«xFWHM))2 )  for  ^  ^  (17) 

where  [H]max  is  the  maximum  in  the  MSISE  density  profile  located  at  an  altitude  altmax,  Dn 
controls  the  diurnal  variability,  and  FWHM  is  the  full  width  at  half  max  of  the  distribution.  For 
Hydrogen,  the  FWHM  employed  is  19.5  km  and  Dn  varies  from  1.0  to  0.5  (day  to  night).  The 
value  of  FWHM  and  its  diurnal  variability  is  chosen  to  best  reproduce  experimental  data.  The 
supplemental  density  is  used  for  altitudes  up  to  the  H  peak  density  near  85  km.  The  resulting  H 
atom  profiles  for  MSISE-90  are  compared  to  the  unmodified  profile  in  Figure  15. 


Figure  15:  A  comparison  of  the  H  atom  profiles  for  two  versions  of  SAG  for  mid  latitude 
summer.  The  SAG  2.0  profiles  are  obtained  using  the  H  corrections. 
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3.1.7  CO  and  CH4 

The  CO  and  CH4  profiles  are  taken  from  the  NRL  database(3)  for  altitudes  up  to  120  km. 
Above  120  km  the  mixing  ratios  are  set  equal  to  the  120  km  values;  this  should  result  in  an  upper 
limit  to  the  expected  ER  radiances,  which  are  extremely  small. 


Figure  16:  The  SAG  water  vapor  profile.  The  four  regions  of  the  profile  are  numbered. 


3.1.8  HzO 

Updates  to  the  H20  profile  based  on  the  URAP(4)  database  may  be  included  from  the 
Parameter  Menu  (see  Table  1).  When  the  URAP  database  is  not  used  then  the  H20  profiles  from 
SAG  version  0.92  are  obtained. 

SAG  0.92  H20  Profiles 

Due  to  problems  in  the  current  NRL  database(3)  for  water  vapor,  an  alternative  empirical 
model  has  been  developed  for  SAG  which  combines  the  models  in  Ref.  1 1  with  the 
recommendation  in  Ref.  42  for  the  upper  troposphere.  Between  0  and  8  km  the  profile  is  a  linear 
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fit  to  the  models  for  relative  humidity,  which  indicate  an  average  humidity  of  approximately 
75%  at  the  ground  and  30%  at  8  km  (Figure  16,  region  number  1).  As  recommended  in  Ref.  42, 
above  8  km  the  mixing  ratio  is  logarithmically  interpolated  to  the  value  of  4  ppmv  (2.5  ppmm)  at 
the  tropopause,  which  we  take  as  the  altitude  of  the  MSISE  temperature  local  minimum  (Figure 
16,  region  number  2).  In  the  unlikely  event  that  the  tropopause  is  below  9  km,  it  is  assumed  to 
be  at  9  km  for  the  purposes  of  this  calculation.  At  40  km  and  above,  the  standard  profile  from 
Ref.  1 1  is  used  (Figure  16,  region  number  4).  Between  the  tropopause  (4  ppm  mixing  ratio)  and 
40  km  (approximately  5  ppm  mixing  ratio)  the  mixing  ratio  is  logarithmically  interpolated 
(Figure  16,  region  number  3). 


Latitude 

Figure  17:  Water  vapor  number  density  differences  between  SAG  0.92  and  SAG  2.0  using 
NRLMSISE-00  and  URAP  options  for  spring,  mid-latitude,  and  moderate-activity 
conditions.  The  figure  shows  the  relative  difference:  ([E^Oko  -  [H2OJ0.92VtH2OJ0.92, 
where  [H2O]  is  the  water  number  density  and  the  subscript  indicates  the  SAG  version. 


SAG  2.0  H20  Profiles 

Updates  to  the  H20  profile  based  on  the  URAP(4)  water  vapor  climatology  may  be  included 
from  the  parameter  menu  (see  Table  1).  The  altitude  range  of  the  URAP  model  is  from  about  16 
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km  to  85  km.  For  higher  altitudes  the  standard  profile  from  D.  Grantham  et  alm  is  used.  The 
empirical  models  employed  by  SAG  for  altitudes  between  0  km  and  the  tropopause  have  been 
retained  (Figure  16,  numbers  1  and  2).  If  the  URAP  water  vapor  climatology  is  not  selected  then 
the  original  SAG  0.92  water  density  profiles  are  obtained  when  the  MSISE-90  option  is  selected. 
A  comparison  of  the  SAG  0.92  water  vapor  density  and  the  updated  SAG  2.0  water  vapor 
density  using  the  URAP  climatology  is  shown  in  Figure  17. 

3.2  Additional  Species  for  SAMM 

The  species  N20,  HN03,  and  N02  are  taken  from  the  NRL  database.(3)  These  profiles 
include  seasonal  and  latitude  variabilities,  but  do  not  address  the  diurnal  variability,  which  in  the 
case  of  N02  is  very  large  due  to  daytime  depletion  by  the  N02  +  O  reaction  and  by  photolysis. 
The  species  S02  and  NH3  are  not  covered  by  the  NRL  database(3),  so  the  standard  AFGL 
constituent  profiles  are  used. 

3.3  Temperature 

The  empirical  MSISE  model  developed  by  Hedin(I,2)  provides  a  smooth  extension  of  the 
widely-used  MSISE  upper  atmospheric  model  to  the  lower  atmosphere.  The  accuracy  of  the 
temperature  profiles  above  120  km,  which  are  based  on  mass  spectrometer,  incoherent  scatter 
and  satellite  drag  data,  has  been  discussed  in  Hedin’s  documentation  for  previous  MSISE 
versions.  The  temperature  profiles  in  MSISE-90  have  been  designed  to  follow  standard  lower 
atmosphere  climatology  compilations  quite  closely  (to  within  several  degrees)  under  most 
conditions.  A  comparison  of  the  MSISE  and  NRL  climatology  temperatures  is  provided  in  the 
optional  SAG  output  file  “dumpfildat”  (see  Appendix  B).  In  most  cases,  close  agreement  is 
found.  The  MSISE  temperature  profiles  are  probably  suitable  for  most  SHARC/SAMM 
applications.  It  should  be  noted,  however,  that  ER  radiances  can  be  extremely  temperature- 
sensitive,  particularly  at  low  temperatures.  One  altitude  region  where  the  temperature 
uncertainty  may  be  significant  is  directly  above  the  mesopause,  where  there  is  a  steep  gradient. 
A  comparison  of  the  calculated  temperatures  for  two  MSISE  models  is  shown  in  Figure  18. 
From  the  figure,  it  is  evident  that  there  are  significant  differences  between  both  the  day  and  night 
temperature  profiles. 
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Figure  18:  The  difference  in  temperatures  between  the  MSISE-90  and  NRLMSIS-00  models  for 
spring,  mid-latitude,  and  moderate-activity  conditions.  The  figure  shows  the  absolute 
difference:  (TNrlmsise-oo  -  Tmsise-9o),  where  T  is  the  temperature  and  the  subscript 
indicates  the  MSISE  version. 
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4.  VERIFICATION  AND  VALIDATION 

SAG  atmospheric  profiles  have  been  compared  to  many  other  profiles  from  independent 
sources  over  a  broad  range  of  conditions.  A  sample  of  these  comparisons  for  all  the  chemical 
species  available  in  SAG  is  presented  here  in  Figure  19-27.  In  all  of  the  figures,  SAG  2.0 
atmospheric  profiles  are  represented  by  the  solid  black  lines.  Comparison  of  the  profiles  shows 
good  agreement  between  the  SAG  and  other  profiles. 
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Figure  19:  A  comparison  with  SAG  daytime,  mid-latitude,  springtime  O  atom  profiles. 
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20:  A  comparison  with  SAG  nighttime,  mid-latitude,  springtime  H  atom  profiles. 
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Figure  21 :  A  comparison  with  SAG  daytime,  mid-latitude,  springtime  ozone  profiles. 
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Figure  22:  A  comparison  with  SAG  mid-latitude,  springtime  CO  profiles. 
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Figure  23:  A  comparison  with  SAG  mid-latitude,  springtime  C02  profiles. 
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Figure  24:  A  comparison  with  SAG  mid-latitude,  springtime  H2O  profiles. 


Figure  25:  A  comparison  with  SAG  mid-latitude,  springtime  OH  profiles. 
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Figure  26:  A  comparison  of  SAG  and  AFGL  atmospheric  profiles  for  mid-latitude  summer. 
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Figure  27:  A  comparison  of  SAG  and  U.S.  Standard  Atmospheric  profiles 


5.  CONCLUSIONS 


The  SHARC/SAMM  Atmosphere  Generator  utilizes  a  combination  of  empirical  models  to 
generate  atmospheric  profiles  for  IR  radiance  calculations  that  properly  account  for  solar 
terminator  and  other  systematic  variability  effects.  The  goal  in  developing  SAG  has  been  to 
provide  reasonably  accurate  temperature  and  species  profiles  with  an  easy-to-use,  rapid¬ 
executing  program.  It  is  anticipated  that  SAG  will  evolve  with  the  radiation  codes  to  keep  up 
with  new  observations,  theoretical  advances,  and  applications  requirements.  Comparisons  of 
available  field  data  with  radiance  predictions  will  comprise  an  important  parallel  effort  and  will 
provide  valuable  tests  of  the  radiation  codes  and  the  atmospheric  profile  models,  contributing  to 
our  overall  understanding  of  atmospheric  processes. 
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APPENDIX  A:  SUMMARY  OF  CHANGES 
A.l  Summary  of  Changes  from  SAG  0.92  to  SAG  2.0 

The  SAG  2.0  program  now: 

•  Includes  the  patched  MSISE90  database  release,  which  corrects  an  error  in 

SUBROUTINE  GHP6  when  using  METER6(.TRUE.).  This  error  does  not  affect  SAG 
output  but  was  corrected  to  aid  future  code  development. 

•  Includes  the  NRLMSISE-00  database. 

•  Offers  a  choice  of 

■  using  either  MSISE-90  or  NRLMSISE-00  database. 

■  including  SNOE  database  for  NO  density. 

■  including  URAP  H2O  database. 

■  using  revised  O  atom  profiles. 

■  using  revised  H  atom  profiles. 

■  using  either  Geodetic  or  Geocentric  Latitude  input. 

•  Is  a  FORTRAN  SUBROUTINE. 

•  Generates  a  default  set  of  altitudes  if  an  altitude  file  is  not  found. 

•  Executes  when  the  input  parameter  file  is  not  found. 

•  Includes  a  quality  assurance  routine  to  test  SAG  installations. 

•  Includes  a  new  atmosphere  output  file  format  for  the  file  ATMOUT.DAT. 

•  Has  a  centralized  error  handling  routine. 
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APPENDIX  B:  PROGRAM  AND  I/O  STRUCTURE 


B.l  Main  Program 

The  main  subroutine,  named  SAG,  serves  to  call  two  subroutines,  the  menu  subroutine  for 
gathering  input  information  and  the  atmosphere  generator  subroutine  which  uses  this 
information.  These  are  described  below. 

B.2  Menu  Subroutine 

The  menu  subroutine,  named  SETDRI,  gathers  information  on  the  day  of  the  year,  universal 
time  or  local  solar  time,  latitude,  longitude,  F10.7  daily  solar  flux  for  previous  day  and  F10.7A 
monthly  average  solar  flux,  daily  Ap,  the  atmosphere  file  name,  and  the  comment  line  appearing 
in  the  atmosphere  file.  At  the  beginning  of  execution,  information  from  the  last  run  is  entered 
from  file  “drivstor.dat”.  At  the  end  of  execution,  the  “drivstor.dat”  data  are  updated. 

B.3  Atmosphere  Generator  Subroutine 

The  atmosphere  generator  subroutine,  named  DRI,  creates  the  atmosphere  file.  Input  files 
opened  are  the  appropriate  “drivht”  file,  which  contains  the  altitude  list  and  the  SHARC/SAMM 
switch,  and  the  NRL  database(3)  files  (see  below).  DRI  calls  on  the  MSISE-90/NRLMSISE-00 
subroutine  GTD6/GTD7. 

In  addition  to  creating  the  SHARC/SAMM  atmosphere  file,  the  delivered  version  of  DRI  may 
also  create  additional  output  files  that  are  useful  for  research  purposes  (see  Section  B.5). 

B.4  List  of  Input  Files 

4  input  files:  drivstor.dat  (updated  by  SAG),  drivht.mod,  drivht.sam,  drivht.sha 

10  NRL  data  files  in  subdirectory  NRLDAT:  CH4.NRL,  CO.NRL,  H20.NRL, 

HN03.NRL,  N20.NRL,  N02.NRL,  03D.NRL,  03R.NRL,  OAT.NRL,  TEMP.NRL. 
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1  URAP  data  file  in  subdirectory  UARSDAT:  H20.DAT 
1  SNOE  data  file  in  subdirectory  SNOE:  SNOE.DAT 
Total  files  =  16. 

B.5  Output  Files 

An  example  of  an  atmosphere  output  file  (which  has  a  user-specified  name)  is  given  in 
Appendix  C.  The  SHARC  atmosphere  files  are  identical  to  the  SAMM  atmosphere  files  except 
for  the  omission  of  the  species  N2O,  HNO3,  NO2,  SO2,  and  NH3. 

The  delivered  code  may  also  create  one  or  more  optional  output  files  useful  for  research 
purposes.  These  are: 

(1)  “tape5.dat”,  an  atmosphere  file  in  the  MODTRAN  input  format.  This  file  is  written  only 
if  the  MODTRAN  or  SAMM  species  (H20,C02,03,N20,C0,CH4,02,N0,S02,N02,NH3,HN03) 
are  requested; 

(2)  “norad.dat”,  a  file  containing  an  altitude  profile  of  the  approximate  NO  5.3  pm  volumetric 
emission,  valid  above  100  km; 

(3)  “dumpfil.dat”,  a  file  containing  a  number  of  outputs  including  pressures,  temperatures, 
number  densities,  concentrations,  and  altitudes  listed  in  columns.  Additional  values  from  the 
NRL  database(3)  (temperature  and  high-altitude  ozone),  which  are  not  used  in  the  output 
atmosphere  are  provided  for  comparison  purposes,  “dumpfil.dat”  differs  from  the  other  output 
files  in  that  each  time  the  “Generate  atmosphere”  command  is  given  within  the  SAG  session  the 
“dumpfil.dat”  output  is  appended  rather  than  overwritten. 

(4)  “ATMOUT.DAT”,  a  file  of  the  calculated  atmosphere  in  an  easily  readable  format. 

(5)  “QAERROR.log”,  a  file  generated  if  a  QA  Test  fails,  which  shows  the  failed  comparison. 
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APPENDIX  C:  SAMPLE  INPUTS  AND  OUTPUTS 


\ 


C.l  Sample  Interactive  Session 

This  section  features  an  illustrative  interactive  session  with  SAG  in  which  an  atmosphere  is 
generated  for  SAMM.  The  user  responses  are  contained  in  braces,  {  }. 

In  this  session,  a  default  value  was  entered  for  the  solar  flux  (moderate)  and  the  local  time 
was  specified  as  10  AM.  The  day  of  the  year  remained  as  day  135.  Note  that  after  the  local  time 
was  changed,  updated  values  for  the  solar  angle  and  the  universal  time  appeared  in  the  Parameter 
Menu  listing. 
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SHARC/SAMM  ATMOSPHERE  GENERATOR 
Version  2.0 

MAIN  MENU 

1  Change  header  line  comment  (64  chars  max)  from 
HIGH  LATITUDE  -  WINTER 

2  Change  output  file  name.  Current  name  is 
SAGTEST.ATM 

3  Change  or  inspect  parameter  values 

4  Enter  defaults 

5  Generate  atmosphere 

6  Run  QA  Tests 

7  Quit 

Enter  item  no: 

{4} 

DEFAULT  MENU 
Current  defaults  are 

1  Day  (d)  or  night  (n): 

2  Equatorial  (e),  mid  (m)  or  high  (h)  latitude: 

3  NH  Winter  (w),  spring  (p),  summer  (s)  or  fall  (f): 

4  Low  (1),  medium  (m)  or  high  (h)  solar  flux:  m 

5  Low  (1),  medium  (m)  or  high  (h)  geomagnetic  activity:  m 

6  SHARC  (s)  SAMM  (a)  MODTRAN  (m)  or  custom  (c)  altitudes:  a 

Enter  [n, letter]  or  [0,0]  to  return  to  main  menu: 

{2,m} 


Enter  [n, letter]  or  [0,0]  to  return  to  main  menu: 

{0} 

1  Day  (d)  or  night  (n): 

2  Equatorial  (e),  mid  (m)  or  high  (h)  latitude:  m 

3  NH  Winter  (w),  spring  (p),  summer  (s)  or  fall  (f): 

4  Low  (1),  medium  (m)  or  high  (h)  solar  flux:  m 

5  Low  (1),  medium  (m)  or  high  (h)  geomagnetic  activity:  m 

6  SHARC  (s)  SAMM  (a)  MODTRAN  (m)  or  custom  (c)  altitudes:  a 


MAIN  MENU 

1  Change  header  line  comment  (64  chars  max)  from 
HIGH  LATITUDE  -  WINTER 

2  Change  output  file  name.  Current  name  is 
SAGTEST.ATM 

3  Change  or  inspect  parameter  values 

4  Enter  defaults 

5  Generate  atmosphere 

6  Run  QA  Tests 

7  Quit 
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Enter  item  no: 

{3} 


PARAMETER  MENU 
Current  parameters  are 

1  Daily  Ap  =  15.0 

2  F10.7  (81  day  centered  average  of  F10.7)  =  120.0 

3  F10.7A  =  120.0 

4  Latitude  Type  =  1  Latitude  (- =S)  =  35.2 
Latitude  Type:  Geocentric(O),  Geodetic(l) 

5  Longitude  =  30.0 

6  UT  (sec)  =  79200. 

7  Solar  Time  (hrs)  =  0.00 

8  Month  =  5 

9  Day  =  15 

10  Year  =  2002 

11  Day  Number  =  135 

12  Solar  Zenith  =  125.82 

13  H,  O,  URAP,  SNOE,  NRLMSISE-00  =TTTTT 
Enter  item  number  [n]  or  [0]  to  return  to  main  menu: 

{7} 

item=  7 

Enter  new  [value]: 

{10} 

Solar  Zenith  Angle  =  64.98 
Enter  item  number  [n]  or  [0]  to  return  to  main  menu: 
{0} 

item=  0 

1  Daily  Ap=  15.0 

2  F10.7  (81  day  centered  average  of  F10.7)  =  120.0 

3  F10.7A  =  120.0 

4  Latitude  Type  =  1  Latitude  (-  =S)  =  35.2 
Latitude  Type:  Geocentric(O),  Geodetic(l) 

5  Longitude  =  30.0 

6  UT  (sec)  =  28800. 

7  Solar  Time  (hrs)  =  10.00 

8  Month  =  5 

9  Day  =  15 

10  Year  =  2002 

1 1  Day  Number  =  135 

12  Solar  Zenith  =  30.45 

13  H,  O,  URAP,  SNOE,  NRLMSISE-00  =  TTTTT 
MAIN  MENU 

1  Change  header  line  comment  (64  chars  max)  from 
HIGH  LATITUDE  -  WINTER 

2  Change  output  file  name.  Current  name  is 
SAGTEST.ATM 

3  Change  or  inspect  parameter  values 

4  Enter  defaults 

5  Generate  atmosphere 


6  Run  QA  Tests 

7  Quit 


Enter  item  no: 

<5> 

sunrise  LST,  sunset  LST=  7.140  16.860 

SHARC/SAMM  ATMOSPHERE  GENERATOR 
Version  2.0 

MAIN  MENU 

1  Change  header  line  comment  (64  chars  max)  from 
HIGH  LATITUDE  -  WINTER 

2  Change  output  file  name.  Current  name  is 
SAGTEST.ATM 

3  Change  or  inspect  parameter  values 

4  Enter  defaults 

5  Generate  atmosphere 

6  Run  QA  Tests 

7  Quit 

Enter  item  no: 

{7} 


C.2  Atmosphere  Output  File  (for  SAMM)  from  the  Interactive  Session 

ATMOSPHERE  NAME  HIGH  LATITUDE  -  WINTER 

SAGTEST.ATM 

END 

NUMBER  OF  LAYERS 

103 

END 

DAY-NIGHT  VARIABLE  AND  EXOATMOSPHERIC  TEMPERATURE 
DAY  1017.81 
END 
SPECIES 

02  O  C02  N2  NO  03  H  OH  CO  H20  CH4  HN03  N20  NH3  N02  S02 
END 

ALTITUDES 


0.00 

1.00 

2.00  ■ 

1.00  4 

.00 

5.00 

6.00 

7.00  8.00  9 

.00 

10.00 

11.00 

12.00 

13.00 

14.00 

15.00 

16.00 

17.00 

18.00 

19.00 

20.00 

21.00 

22.00 

23.00 

24.00 

26.00 

28.00 

30.00 

32.00 

34.00 

36.00 

38.00 

40.00 

42.00 

44.00 

46.00 

48.00 

50.00 

52.00 

54.00 

56.00 

58.00 

60.00 

62.00 

64.00 

66.00 

68.00 

70.00 

72.00 

74.00 

76.00 

78.00 

80.00 

82.00 

84.00 

86.00 

88.00 

90.00 

92.00 

94.00 

96.00 

98.00 

100.00 

102.00  104.00 

106.00  108.00  110.00  112.00  114.00 
116.00  118.00  120.00  122.00  124.00 
126.00  128.00  130.00  132.00  134.00 
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136.00  138.00  140.00  142.00  144.00 
146.00  148.00  150.00  160.00  170.00 
180.00  190.00  200.00  210.00  220.00 
230.00  240.00  250.00  260.00  270.00 
280.00  290.00  300.00 
END 

TEMPERATURES 

0.2932E+03  0.2890E+03  0.2837E+03  0.2776E+03  0.2709E+03 
0.2639E+03  0.2568E+03  0.2497E+03  0.2427E+03  0.236 1E+03 
0.2299E+03  0.2241E+03  0.2189E+03  0.2146E+03  0.2113E+03 
0.209 1E+03  0.2081E+03  0.2081E+03  0.2089E+03  0.2102E+03 
0.2118E+03  0.2134E+03  0.2151E+03  0.2169E+03  0.2187E+03 
0.2226E+03  0.2267E+03  0.2312E+03  0.236 1E+03  0.2417E+03 
0.248 1E+03  0.2550E+03  0.2617E+03  0.2675E+03  0.2716E+03 
0.2732E+03  0.2720E+03  0.2688E+03  0.2640E+03  0.2584E+03 
0.2525E+03  0.2467E+03  0.2408E+03  0.235 1E+03  0.2296E+03 
0.2241E+03  0.2189E+03  0.2138E+03  0.2090E+03  0.2040E+03 
0.1983E+03  0.1923E+03  0.1865E+03  0.1813E+03  0.1770E+03 
0.1738E+03  0.1718E+03  0.1714E+03  0.1727E+03  0.1757E+03 
0.1806E+03  0.1874E+03  0.1965E+03  0.2079E+03  0.2219E+03 
0.238 1E+03  0.2562E+03  0.2750E+03  0.2934E+03  0.3110E+03 
0.3284E+03  0.3462E+03  0.3658E+03  0.3889E+03  0.4179E+03 
0.4478E+03  0.4762E+03  0.503 1E+03  0.5287E+03  0.5530E+03 
0.5761E+03  0.5980E+03  0.6189E+03  0.6387E+03  0.6574E+03 
0.6753E+03  0.6923E+03  0.7084E+03  0.7776E+03  0.8312E+03 
0.8727E+03  0.9049E+03  0.9299E+03  0.9493E+03  0.9644E+03 
0.976 1E+03  0.9852E+03  0.9923E+03  0.9979E+03  0.1002E+04 
0.1006E+04  0.1008E+04  0.1010E+04 
END 

N2  DENSITIES 

0.1977E+20  0.1784E+20  0.1613E+20  0.1460E+20  0.1321E+20 
0.1 194E+20  0.1077E+20  0.9682E+19  0.8672E+19  0.7734E+19 
0.6865E+19  0.6062E+19  0.5322E+19  0.4641E+19  0.4018E+19 
0.3455E+19  0.2950E+19  0.2505E+19  0.2121E+19  0.1793E+19 
0.1515E+19  0.1282E+19  0.1086E+19  0.9206E+18  0.7815E+18 
0.5650E+18  0.4104E+18  0.2995E+18  0.2196E+18  0.1617E+18 
0.1196E+18  0.8898E+17  0.6680E+17  0.5066E+17  0.3888E+17 
0.3020E+17  0.2370E+17  0.1872E+17  0.1481E+17  0.1170E+17 
0.9209E+16  0.7208E+16  0.5608E+16  0.4336E+16  0.333 1E+16 
0.2541E+16  0.1924E+16  0.1447E+16  0.1080E+16  0.8002E+15 
0.5914E+15  0.4337E+15  0.3148E+15  0.2257E+15  0.1597E+15 
0.1116E+15  0.7704E+14  0.5263E+14  0.3567E+14  0.2408E+14 
0.1625E+14  0.1 101E+14  0.7496E+13  0.5146E+13  0.3574E+13 
0.2518E+13  0.1806E+13  0.1322E+13  0.9895E+12  0.7556E+12 
0.5862E+12  0.4602E+12  0.3642E+12  0.2893E+12  0.2300E+12 
0.1853E+12  0.1519E+12  0.1263E+12  0.1063E+12  0.9038E+11 
0.7756E+11  0.6708E+11  0.5843E+11  0.5121E+11  0.4513E+11 
0.3997E+11  0.3555E+11  0.3175E+11  0.1896E+1 1  0.1201E+1 1 
0.7931E+10  0.5388E+10  0.3738E+10  0.2634E+10  0.1878E+10 
0.1351E+10  0.979 1E+09  0.7134E+09  0.5222E+09  0.3837E+09 
0.2827E+09  0.2089E+09  0.1546E+09 
END 

02  DENSITIES 

0.5303E+19  0.4785E+19  0.4328E+19  0.3917E+19  0.3545E+19 
0.3204E+19  0.2889E+19  0.2597E+19  0.2326E+19  0.2075E+19 
0.1842E+19  0.1626E+19  0.1428E+19  0.1245E+19  0.1078E+19 
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0.9268E+18  0.7913E+18  0.6721E+18  0.5690E+18  0.4810E+18 
0.4065E+18  0.3439E+18  0.2913E+18  0.2470E+18  0.2096E+18 
0.1516E+18  0.1 101E+18  0.8036E+17  0.5892E+17  0.4338E+17 
0.3208E+17  0.2387E+17  0.1792E+17  0.1359E+17  0.1043E+17 
0.8102E+16  0.6359E+16  0.5021E+16  0.3973E+16  0.3139E+16 
0.247 1E+ 16  0.1934E+16  0.1505E+16  0.1 163E+16  0.8908E+15 
0.6768E+15  0.5106E+15  0.3823E+15  0.2842E+15  0.2098E+15 
0.1544E+15  0.1 127E+15  0.8128E+14  0.5786E+14  0.4058E+14 
0.2805E+14  0.1911E+14  0.1284E+14  0.8534E+13  0.5627E+13 
0.3693E+13  0.2421E+13  0.1588E+13  0.1046E+13  0.6943E+12 
0.4660E+12  0.3174E+12  0.2202E+12  0.1560E+12  0.1126E+12 
0.8252E+1 1  0.61 16E+1 1  0.4570E+11  0.3431E+11  0.2581E+11 
0.1973E+1 1  0.1538E+1 1  0.1218E+11  0.9799E+10  0.7985E+10 
0.6584E+10  0.5487E+10  0.4616E+10  0.3918E+10  0.3351E+10 
0.2887E+10  0.2503E+10  0.2183E+10  0.1 182E+10  0.6960E+09 
0.4322E+09  0.278 1E+09  0.1834E+09  0.1231E+09  0.8377E+08 
0.5758E+08  0.3988E+08  0.2780E+08  0.1948E+08  0.1370E+08 
0.9667E+07  0.6842E+07  0.4854E+07 
END 

O  DENSITIES 

0.1308E+03  0.2443E+03  0.4509E+03  0.8229E+03  0.1484E+04 
0.2648E+04  0.4668E+04  0.8137E+04  0.1402E+05  0.2389E+05 
0.4023E+05  0.6699E+05  0.1 103E+06  0.1794E+06  0.2887E+06 
0.4591E+06  0.7219E+06  0.1 122E+07  0.1724E+07  0.2620E+07 
0.3934E+07  0.5841E+07  0.8574E+07  0.1244E+08  0.1785E+08 
0.3550E+08  0.6743E+08  0.1223E+09  0.2120E+09  0.3510E+09 
0.5551E+09  0.8385E+09  0.1210E+10  0.1667E+10  0.2195E+10 
0.2760E+10  0.3316E+10  0.3804E+10  0.4170E+10  0.4365E+10 
0.4380E+10  0.4296E+10  0.4135E+10  0.3903E+10  0.3615E+10 
0.3284E+10  0.2927E+10  0.2559E+10  0.2195E+10  0.1870E+10 
0.1662E+10  0.1949E+10  0.4052E+10  0.1 144E+1 1  0.3030E+1 1 
0.6706E+11  0.1226E+12  0.1881E+12  0.2484E+12  0.2900E+12 
0.3072E+12  0.3026E+12  0.2825E+12  0.2540E+12  0.223 1E+12 
0.1935E+12  0.1671E+12  0.1448E+12  0.1266E+12  0.11 18E+12 
0.9938E+11  0.8883E+1 1  0.7957E+11  0.7118E+11  0.6086E+11 
0.5243E+1 1  0.4570E+1 1  0.4024E+11  0.3573E+11  0.3197E+11 
0.2879E+1 1  0.2607E+1 1  0.2373E+11  0.2170E+11  0.1993E+11 
0.1837E+11  0.1699E+11  0.1576E+11  0.1124E+11  0.8398E+10 
0.6476E+10  0.5106E+10  0.4091E+10  0.3318E+10  0.2715E+10 
0.2237E+10  0.1853E+10  0.1542E+10  0.1287E+10  0.1077E+10 
0.9029E+09  0.7586E+09  0.6383E+09 
END 

C02  DENSITIES 

0.8607E+16  0.7767E+16  0.7024E+16  0.6358E+16  0.5754E+16 
0.5200E+16  0.4689E+16  0.4216E+16  0.3776E+16  0.3368E+16 
0.2989E+16  0.2640E+16  0.2317E+16  0.2021E+16  0.1750E+16 
0.1504E+16  0.1284E+16  0.1091E+16  0.9236E+15  0.7807E+15 
0.6598E+15  0.5582E+15  0.4728E+15  0.4008E+15  0.3403E+15 
0.2460E+15  0.1787E+15  0.1304E+15  0.9563E+14  0.7042E+14 
0.5206E+14  0.3875E+14  0.2908E+14  0.2206E+14  0.1693E+14 
0.1315E+14  0.1032E+14  0.8149E+13  0.6448E+13  0.5095E+13 
0.4010E+13  0.3139E+13  0.2442E+13  0.1888E+13  0.1450E+13 
0.1 105E+13  0.8369E+12  0.6290E+12  0.4692E+12  0.3475E+12 
0.2565E+12  0.1879E+12  0.1361E+12  0.9730E+1 1  0.6858E+1 1 
0.4763 E+ 11  0.3261E+11  0.2201E+11  0.1384E+11  0.8614E+10 
0.5324E+10  0.3274E+10  0.2007E+10  0.1229E+10  0.7551E+09 


-48- 


0.4674E+09  0.2930E+09  0.1869E+09  0.1218E+09  0.8095E+08 
0.5466E+08  0.3737E+08  0.2577E+08  0.1874E+08  0.1365E+08 
0.1012E+08  0.7667E+07  0.5914E+07  0.4634E+07  0.3681E+07 
0.2959E+07  0.2403E+07  0.1970E+07  0.1628E+07  0.1356E+07 
0.1136E+07  0.9580E+06  0.8119E+06  0.3829E+06  0.1953E+06 
0.1050E+06  0.5864E+05  0.3365E+05  0.1971E+05  0.1172E+05 
0.7054E+04  0.4286E+04  0.2624E+04  0.1616E+04  0.9999E+03 
0.6212E+03  0.3873E+03  0.2422E+03 
END 

H  DENSITIES 

0.1000E-04  0.1000E-04  0.1000E-04  0.1000E-04  0.1000E-04 
0.1000E-04  0. 1 000E-04  0.1000E-04  0.1000E-04  0.1000E-04 
0.1000E-04  0.1000E-04  0.1000E-04  0.1000E-04  0.1000E-04 
0.1000E-04  0.1000E-04  0.1000E-04  0.1000E-04  0.2330E-04 
0.5697E-04  0.1374E-03  0.3268E-03  0.7665E-03  0.1773E-02 
0.9103E-02  0.442 IE-01  0.2031E+00  0.8830E+00  0.3631E+01 
0.1413E+02  0.5200E+02  0.1811E+03  0.5966E+03  0.1859E+04 
0.5482E+04  0.1529E+05  0.4036E+05  0.1008E+06  0.2380E+06 
0.5318E+06  0.1 124E+07  0.2249E+07  0.4254E+07  0.7616E+07 
0.1290E+08  0.2066E+08  0.3132E+08  0.4491E+08  0.6093E+08 
0.7819E+08  0.9494E+08  0.1091E+09  0.1 185E+09  0.1260E+09 
0.1 132E+09  0.8872E+08  0.6665E+08  0.4991E+08  0.3794E+08 
0.2950E+08  0.2353E+08  0.1922E+08  0.1605E+08  0.1361E+08 
0.1166E+08  0.1005E+08  0.8672E+07  0.7457E+07  0.6387E+07 
0.5449E+07  0.4633E+07  0.3935E+07  0.3348E+07  0.2871E+07 
0.248 1E+07  0.2157E+07  0.1886E+07  0.1657E+07  0.1463E+07 
0.1297E+07  0.1156E+07  0.1034E+07  0.9284E+06  0.8375E+06 
0.7588E+06  0.6904E+06  0.6309E+06  0.4279E+06  0.3197E+06 
0.2588E+06  0.2228E+06  0.2006E+06  0.1863E+06  0.1766E+06 
0.1699E+06  0.1649E+06  0.1610E+06  0.1578E+06  0.1551E+06 
0.1527E+06  0.1505E+06  0.1485E+06 
END 

NO  DENSITIES 

0.7592E+10  0.6851E+10  0.6196E+10  0.5608E+10  0.5075E+10 
0.4587E+10  0.4136E+10  0.3718E+10  0.3331E+10  0.2970E+10 
0.2637E+10  0.2328E+10  0.2044E+10  0.1776E+10  0.1518E+10 
0.1252E+10  0.1012E+10  0.8083E+09  0.6517E+09  0.5601E+09 
0.4947E+09  0.4546E+09  0.4267E+09  0.4243E+09  0.4512E+09 
0.6363E+09  0.7659E+09  0.9395E+09  0.1 126E+10  0.1232E+10 
0.1217E+10  0.1103E+10  0.9578E+09  0.7625E+09  0.5863E+09 
0.4414E+09  0.3294E+09  0.2468E+09  0.1937E+09  0.1519E+09 
0.1191E+09  0.9320E+08  0.7252E+08  0.5651E+08  0.4372E+08 
0.3420E+08  0.2705E+08  0.2123E+08  0.1812E+08  0.1535E+08 
0.1342E+08  0.1206E+08  0.1072E+08  0.121 1E+08  0.1368E+08 
0.1545E+08  0.1746E+08  0.1972E+08  0.2228E+08  0.2518E+08 
0.2844E+08  0.3213E+08  0.3630E+08  0.8448E+08  0.1052E+09 
0.1237E+09  0.1352E+09  0.1370E+09  0.1297E+09  0.1165E+09 
0.1012E+09  0.8664E+08  0.7404E+08  0.6368E+08  0.5573E+08 
0.4999E+08  0.4592E+08  0.4292E+08  0.3446E+08  0.2902E+08 
0.247 1E+08  0.2097E+08  0.1780E+08  0.1525E+08  0.1307E+08 
0.1116E+08  0.9535E+07  0.8178E+07  0.5601E+07  0.3834E+07 
0.2622E+07  0.1781E+07  0.1210E+07  0.8204E+06  0.5560E+06 
0.3763E+06  0.2545E+06  0.1718E+06  0.1 159E+06  0.7806E+05 
0.5250E+05  0.3525E+05  0.2362E+05 
END 

03  DENSITIES 
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0.1003E+13  0.1066E+13  0.1137E+13  0.1063E+13  0.9942E+12 
0.9314E+12  0.8707E+12  0.8371E+12  0.8018E+12  0.8730E+12 
0.9460E+ 12  0.11 39E+ 13  0.1 364E+ 1 3  0. 1 628E+ 1 3  0. 1 929E+ 1 3 
0.2256E+13  0.2621E+13  0.2854E+13  0.3098E+13  0.3356E+13 
0.3635E+13  0.3942E+13  0.4281E+13  0.4654E+13  0.5065E+13 
0.4284E+13  0.3577E+13  0.3001E+13  0.2269E+13  0.1723E+13 
0.1246E+13  0.8597E+12  0.5985E+12  0.3776E+12  0.2410E+12 
0.1555E+12  0.1012E+12  0.6630E+1 1  0.4477E+11  0.3019E+11 
0.1990E+1 1  0.1280E+1 1  0.8189E+10  0.4993E+10  0.3023E+10 
0.1775E+10  0.101 1E+10  0.5720E+09  0.3716E+09  0.2397E+09 
0.1712E+09  0.1 123E+09  0.5757E+08  0.6015E+08  0.6463E+08 
0.6495E+08  0.5959E+08  0.4366E+08  0.2634E+08  0.1364E+08 
0.6333E+07  0.2752E+07  0.1164E+07  0.4976E+06  0.2224E+06 
0.1059E+06  0.5394E+05  0.2928E+05  0.1686E+05  0.1015E+05 
0.6307E+04  0.4000E+04  0.2559E+04  0.1629E+04  0.9892E+03 
0.6164E+03  0.3999E+03  0.2687E+03  0.1862E+03  0.1325E+03 
0.9665E+02  0.7199E+02  0.5463E+02  0.4216E+02  0.3302E+02 
0.262 1E+02  0.2105E+02  0.1709E+02  0.6823E+01  0.3148E+01 
0.1587E+01  0.8471E+00  0.4695E+00  0.2669E+00  0.1545E+00 
0.9053E-01  0.5353E-01  0.3186E-01  0.1906E-01  0.1 144E-01 
0.6884E-02  0.4152E-02  0.2509E-02 
END 

OH  DENSITIES 

0.1113E+07  0.1005E+07  0.9087E+06  0.8225E+06  0.7443E+06 
0.6727E+06  0.6066E+06  0.5466E+06  0.4940E+06  0.4515E+06 
0.41 13E+06  0.3725E+06  0.3366E+06  0.3083E+06  0.2907E+06 
0.2985E+06  0.31 15E+06  0.3336E+06  0.3530E+06  0.3765E+06 
0.4190E+06  0.5580E+06  0.7075E+06  0.8945E+06  0.1161E+07 
0.2198E+07  0.31 14E+07  0.4487E+07  0.7723E+07  0.1207E+08 
0.1699E+08  0.2147E+08  0.2437E+08  0.2453E+08  0.2320E+08 
0.2077E+08  0.1812E+08  0.1617E+08  0.1469E+08  0.1333E+08 
0.1270E+08  0.1271E+08  0.1264E+08  0.1321E+08  0.1371E+08 
0.1384E+08  0.1361E+08  0.1328E+08  0.1182E+08  0.6944E+07 
0.4080E+07  0.2397E+07  0.1408E+07  0.5358E+06  0.2200E+06 
0.8637E+05  0.3317E+05  0.1183E+05  0.4115E+04  0.1439E+04 
0.5187E+03  0.1963E+03  0.7944 E+02  0.3490E+02  0.1681E+02 
0.8847E+01  0.5006E+01  0.2981E+01  0.1833E+01  0.1149E+01 
0.7285E+00  0.4654E+00  0.2985E+00  0.1918E+00  0.1246E+00 
0.8252E-01  0.5604E-01  0.3895E-01  0.2766E-01  0.2003E-01 
0.1476E-01  0.1 107E-01  0.8421E-02  0.6500E-02  0.5083E-02 
0.4024E-02  0.3222E-02  0.2607E-02  0.1036E-02  0.4927E-03 
0.2663E-03  0.1575E-03  0.9914E-04  0.6506E-04  0.4389E-04 
0.3016E-04  0.2097E-04  0.1469E-04  0.1034E-04  0.7299E-05 
0.5164E-05  0.3658E-05  0.2593E-05 
END 

CO  DENSITIES 

0.2784E+13  0.2467E+13  0.2210E+13  0.1973E+13  0.1760E+13 
0.1575E+13  0.1407E+13  0.1253E+13  0.1117E+13  0.9860E+12 
0.8667E+12  0.7574E+12  0.6514E+12  0.5383E+12  0.4380E+12 
0.3236E+12  0.2310E+12  0.1546E+12  0.9559E+11  0.7383E+11 
0.5658E+1 1  0.4368E+1 1  0.3417E+11  0.2601E+11  0.2004E+11 
0.1247F.+1 1  0.7682E+10  0.4602E+10  0.4034E+10  0.3456E+10 
0.29 1 9E+ 1 0  0.2446E+ 1 0  0.204 1 E+ 1 0  0. 1 844E+ 1 0  0. 1 64 1 E+ 1 0 
0.1640E+10  0.1721E+10  0.1701E+10  0.1889E+10  0.1921E+10 
0.1991E+10  0.2043E+10  0.1967E+10  0.1978E+10  0.1870E+10 
0.1742E+10  0.1596E+10  0.1407E+10  0.1362E+10  0.1240E+10 
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0.1 194E+10  0.1 158E+10  0.1045E+10  0.1007E+10  0.8951E+09 
0.7708E+09  0.6449E+09  0.5175E+09  0.4202E+09  0.3307E+09 
0.2591E+09  0.2027E+09  0.1567E+09  0.1212E+09  0.9375E+08 
0.7242E+08  0.5625E+08  0.4441E+08  0.3519E+08  0.2842E+08 
0.2313E+08  0.1891E+08  0.1560E+08  0.1258E+08  0.1008E+08 
0.8198E+07  0.6779E+07  0.5688E+07  0.4831E+07  0.4147E+07 
0.3593E+07  0.3138E+07  0.2760E+07  0.2443E+07  0.2175E+07 
0.1947E+07  0.1750E+07  0.1580E+07  0.9973E+06  0.6707E+06 
0.4715E+06  0.3422E+06  0.2545E+06  0.1930E+06  0.1485E+06 
0.1158E+06  0.9121E+05  0.7249E+05  0.5806E+05  0.4681E+05 
0.3796E+05  0.3094E+05  0.2533E+05 
END 

H20  DENSITIES 

0.4345E+18  0.3130E+18  0.2073E+18  0.1271E+18  0.7269E+17 
0.3902E+17  0.1979E+17  0.9531E+16  0.4373E+16  0.2197E+16 
0.1099E+16  0.5467E+15  0.2704E+15  0.1329E+15  0.6482E+14 
0.3140E+14  0.1510E+14  0.1236E+14  0.1039E+14  0.9051E+13 
0.8063E+13  0.7143E+13  0.6160E+13  0.5212E+13  0.4384E+13 
0.3167E+13  0.2416E+13  0.1864E+13  0.1414E+13  0.1061E+13 
0.8017E+12  0.6108E+12  0.4666E-rl2  0.3580E+12  0.2792E+12 
0.2202E+12  0.1747E+12  0.1406E+12  0.1 130E+12  0.8921E+1 1 
0.6942E+11  0.5277E+11  0.3984E+11  0.3079E+11  0.2356E+11 
0.1810E+11  0.1344E+11  0.9018E+10  0.5787E+10  0.3652E+10 
0.2205E+10  0.1270E+10  0.7876E+09  0.4794E+09  0.2890E+09 
0.1715E+09  0.9879E+08  0.5632E+08  0.3180E+08  0.1789E+08 
0.1039E+08  0.6243E+07  0.3776E+07  0.2434E+07  0.1589E+07 
0.1046E+07  0.6972E+06  0.4749E+06  0.3366E+06  0.2437E+06 
0.1785E+06  0.1318E+06  0.9824E+05  0.7922E+05  0.635 1E+05 
0.5163E+05  0.4270E+05  0.3583E+05  0.3043E+05  0.2612E+05 
0.2263E+05  0.1976E+05  0.1739E+05  0.1539E+05  0.1370E+05 
0.1226E+05  0.1102E+05  0.9952E+04  0.6282E+04  0.4225E+04 
0.2970E+04  0.2155E+04  0.1603E+04  0.1215E+04  0.9356E+03 
0.7293E+03  0.5745E+03  0.4566E+03  0.3657E+03  0.2948E+03 
0.2391E+03  0.1949E+03  0.1596E+03 
END 

CH4  DENSITIES 

0.4366E+14  0.3888E+14  0.3470E+14  0.3097E+14  0.2765E+14 
0.2464E+14  0.2191E+14  0.1941E+14  0.1714E+14  0.1506E+14 
0.1317E+14  0.1146E+14  0.9901E+13  0.8500E+13  0.7241E+13 
0.6125E+13  0.5144E+13  0.4296E+13  0.3576E+13  0.2972E+13 
0.2468E+13  0.2050E+13  0.1705E+13  0.1419E+13  0.1182E+13 
0.8215E+12  0.5730E+12  0.4009E+12  0.2671E+12  0.1769E+12 
0.1 166E+12  0.7655E+1 1  0.4980E+11  0.3473E+11  0.2431E+11 
0.1747E+11  0.1292E+1 1  0.9571E+10  0.7153E+10  0.5320E+10 
0.3943E+10  0.2909E+10  0.2126E+10  0.1551E+10  0.1120E+10 
0.8035E+09  0.5734E+09  0.4047E+09  0.2848E+09  0.1983E+09 
0.1378E+09  0.9518E+08  0.6478E+08  0.4386E+08  0.292 1E+08 
0.1919E+08  0.1244E+08  0.7954E+07  0.5088E+07  0.3233E+07 
0.2056E+07  0.1313E+07  0.8404E+06  0.5468E+06  0.3589E+06 
0.2390E+06  0. 1 620E+06  0.1 1 18E+06  0.7939E+05  0.573 7E+05 
0.4237E+05  0.3191E+05  0.2419E+05  0.1950E+05  0.1564E+05 
0.1271E+05  0.1051E+05  0.8820E+04  0.7491E+04  0.6430E+04 
0.5571E+04  0.4865E+04  0.4280E+04  0.3789E+04  0.3373E+04 
0.3019E+04  0.2714E+04  0.2450E+04  0.1546E+04  0.1040E+04 
0.7310E+03  0.5306E+03  0.3947E+03  0.2992E+03  0.2303E+03 
0.1795E+03  0.1414E+03  0.1124E+03  0.9003E+02  0.7258E+02 
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0.5886E+02  0.4797E+02  0.3928E+02 
END 

HN03  DENSITIES 

0.1204E+11  0.1202E+11  0.1205E+11  0.1207E+11  0.1210E+11 
0.1211E+11  0.1208E+1 1  0.1204E+11  0.1195E+11  0.1178E+11 
0.1 159E+1 1  0.1132E+11  0.1100E+11  0.1063E+11  0.1020E+11 
0.9706E+10  0.9172E+10  0.8632E+10  0.8090E+10  0.7577E+10 
0.7944E+10  0.8037E+10  0.7531E+10  0.6963E+10  0.5966E+10 
0.3880E+10  0.2315E+10  0.1388E+10  0.7799E+09  0.4401E+09 
0.2257E+09  0.1055E+09  0.4974E+08  0.2371E+08  0.1143E+08 
0.5578E+07  0.2750E+07  0.1364E+07  0.6779E+06  0.3366E+06 
0.1664E+06  0.8181E+05  0.3998E+05  0.1941E+05  0.9359E+04 
0.4480E+04  0.2130E+04  0.1005E+04  0.4708E+03  0.2190E+03 
0.1016E+03  0.4675E+02  0.2128E+02  0.9575E+01  0.4251E+01 
0.1862E+01  0.8058E+00  0.3451E+00  0.1467E+00  0.6217E-01 
0.2637E-01  0.1123E-01  0.4810E-02  0.2078E-02  0.9090E-03 
0.4037E-03  0.1827E-03  0.8446E-04  0.4000E-04  0.1935E-04 
0.9526E-05  0.4752E-05  0.2393E-05  0.1930E-05  0.1547E-05 
0.1258E-05  0.1040E-05  0.8727E-06  0.7412E-06  0.6363E-06 
0.5512E-06  0.4814E-06  0.4235E-06  0.3749E-06  0.3338E-06 
0.2987E-06  0.2685E-06  0.2424E-06  0.1530E-06  0.1029E-06 
0.7234E-07  0.5251E-07  0.3905E-07  0.2961E-07  0.2279E-07 
0.1777E-07  0.1399E-07  0.1 1 12E-07  0.8908E-08  0.7182E-08 
0.5824E-08  0.4747E-08  0.3887E-08 
END 

N20  DENSITIES 

0.7695E+13  0.6881E+13  0.6208E+13  0.5619E+13  0.5085E+13 
0.4596E+13  0.4144E+13  0.3726E+13  0.3337E+13  0.2976E+13 
0.2642E+13  0.2333E+13  0.2048E+13  0.1786E+13  0.1546E+13 
0.1320E+13  0.1094E+13  0.9005E+12  0.7390E+12  0.6046E+12 
0.494 1E+ 12  0.4037E+12  0.3298E+12  0.2694E+12  0.2196E+12 
0. 1455E+12  0.9627E+1 1  0.6399E+11  0.3541E+11  0.1967E+11 
0.9838E+10  0.4439E+10  0.2020E+10  0.1202E+10  0.7230E+09 
0.4685E+09  0.3262E+09  0.2286E+09  0.1605E+09  0.1125E+09 
0.7855E+08  0.5453E+08  0.3763E+08  0.2582E+08  0.1758E+08 
0.1 189E+08  0.7981E+07  0.5318E+07  0.3519E+07  0.2312E+07 
0.1514E+07  0.9840E+06  0.6326E+06  0.4017E+06  0.2518E+06 
0.1557E+06  0.9520E+05  0.5761E+05  0.3460E+05  0.2071E+05 
0.1240E+05  0.7456E+04  0.4510E+04  0.2752E+04  0.1699E+04 
0.1066E+04  0.681 1E+03  0.4447E+03  0.2975E+03  0.2032E+03 
0.1413E+03  0.9956E+02  0.7084E+02  0.5712E+02  0.4580E+02 
0.3723E+02  0.3079E+02  0.2583E+02  0.2194E+02  0.1883E+02 
0. 1 632E+02  0.1425E+02  0.1254E+02  0.1 1 10E+02  0.9880E+01 
0.8842E+01  0.7949E+01  0.7176E+01  0.4530E+01  0.3046E+01 
0.2141E+01  0.1554E+01  0.1 156E+01  0.8764E+00  0.6746E+00 
0.5259E+00  0.4143E+00  0.3292E+00  0.2637E+00  0.2126E+00 
0. 1 724E+00  0. 1 405E+00  0.1151  E+00 
END 

N02  DENSITIES 

0.5589E+09  0.6058E+09  0.6575E+09  0.7142E+09  0.7757E+09 
0.8413E+09  0.9104E+09  0.9822E+09  0.1056E+10  0.1131E+10 
0.1205E+10  0.1277E+10  0.1345E+10  0.1408E+10  0.1463E+10 
0.1509E+10  0.1548E+10  0.1578E+10  0.1603E+10  0.1626E+10 
0.1698E+10  0.1723E+10  0.1921E+10  0.2041E+10  0.2351E+10 
0.2526E+10  0.2518E+10  0.2522E+10  0.2128E+10  0.1804E+10 
0.141 1E+10  0.1022E+10  0.7460E+09  0.5193E+09  0.3657E+09 


0.2485E+09  0.1627E+09  0.1071E+09  0.5419E+08  0.2738E+08 
0.1378E+08  0.6901E+07  0.3435E+07  0.1699E+07  0.8343E+06 
0.4068E+06  0.1970E+06  0.9465E+05  0.4513E+05  0.2137E+05 
0.1009E+05  0.4727E+04  0.2191E+04  0.1004E+04  0.4536E+03 
0.2023E+03  0.8918E+02  0.3891E+02  0.1686E+02  0.7284E+01 
0.3150E+01  0.1368E+01  0.5978E+00  0.2645E+00  0.1185E+00 
0.542 IE-01  0.254 IE-01  0.1217E-01  0.6 160E-02  0.31 84E-02 
0.1460E-02  0.5917E-03  0.2421E-03  0.1953E-03  0.1565E-03 
0.1273E-03  0.1052E-03  0.8830E-04  0.7500E-04  0.6438E-04 
0.5577E-04  0.4871E-04  0.4285E-04  0.3793E-04  0.3377E-04 
0.3022E-04  0.2717E-04  0.2453E-04  0.1548E-04  0.1041E-04 
0.7319E-05  0.5313E-05  0.3951E-05  0.2996E-05  0.2306E-05 
0.1798E-05  0.1416E-05  0.1 125E-05  0.9013E-06  0.7266E-06 
0.5892E-06  0.4803E-06  0.3933E-06 
END 

S02  DENSITIES 

0.7592E+10  0.6257E+10  0.4874E+10  0.3552E+10  0.2470E+10 
0.1804E+10  0.1339E+10  0.1029E+10  0.8005E+09  0.6495E+09 
0.5344E+09  0.4494E+09  0.3815E+09  0.3321E+09  0.2901E+09 
0.2543E+09  0.2171E+09  0.1722E+09  0.1298E+09  0.9113E+08 
0.6189E+08  0.4382E+08  0.3169E+08  0.2440E+08  0.1901E+08 
0.1203E+08  0.7870E+07  0.5139E+07  0.3472E+07  0.2442E+07 
0.1806E+07  0.141 1E+07  0.1163E+07  0.1030E+07  0.9491E+06 
0.907 1E+06  0.8838E+06  0.8530E+06  0.7471E+06  0.6535E+06 
0.5538E+06  0.4539E+06  0.3698E+06  0.2850E+06  0.2181E+06 
0.1605E+06  0.1 136E+06  0.7977E+05  0.5022E+05  0.3140E+05 
0. 1 832E+05  0.9922E+04  0.5318E+04  0.2686E+04  0.1339E+04 
0.6535E+03  0.3127E+03  0.1480E+03  0.6925E+02  0.3225E+02 
0.1501E+02  0.6997E+01  0.3279E+01  0.1548E+01  0.7387E+00 
0.3579E+00  0.1765E+00  0.8889E-01  0.4574E-01  0.2401E-01 
0.1281E-01  0.6926E-02  0.3775E-02  0.2980E-02  0.2354E-02 
0.1887E-02  0.1539E-02  0.1274E-02  0.1067E-02  0.9045E-03 
0.7736E-03  0.6671E-03  0.5795E-03  0.5066E-03  0.4455E-03 
0.3938E-03  0.3496E-03  0.3118E-03  0.1849E-03  0.1167E-03 
0.7677E-04  0.5201E-04  0.3599E-04  0.2530E-04  0.1800E-04 
0.1293E-04  0.9348E-05  0.6800E-05  0.4969E-05  0.3645E-05 
0.2682E-05  0.1979E-05  0.1463E-05 
END 

NH3  DENSITIES 

0.1265E+11  0.1142E+11  0.9562E+10  0.7103E+10  0.4872E+10 
0.31 19E+10  0.2013E+10  0.1225E+10  0.7194E+09  0.3733E+09 
0.1784E+09  0.8460E+08  0.4292E+08  0.1854E+08  0.5710E+07 
0.1977E+07  0.7967E+06  0.3528E+06  0.1819E+06  0.9113E+05 
0.4675E+05  0.3151E+05  0.2391E+05  0.1874E+05  0.1441E+05 
0.7980E+04  0.4555E+04  0.2435E+04  0.1 154E+04  0.4994E+03 
0.2126E+03  0.8386E+02  0.2514E+02  0.7212E+01  0.2279E+01 
0.8269E+00  0.3498E+00  0.1708E+00  0.1154E+00  0.7783E-01 
0.5360E-01  0.3764E-01  0.2628E-01  0.1876E-01  0.1330E-01 
0.9447E-02  0.6720E-02  0.4746E-02  0.3367E-02  0.2373E-02 
0.1673E-02  0.1 174E-02  0.8156E-03  0.5615E-03  0.3815E-03 
0.2567E-03  0.1713E-03  0.1130E-03  0.7408E-04  0.4835E-04 
0.3154E-04  0.2062E-04  0.1355E-04  0.9033E-05  0.6087E-05 
0.4158E-05  0.2888E-05  0.2048E-05  0.1488E-05  0.1103E-05 
0.8317E-06  0.6356E-06  0.4899E-06  0.3867E-06  0.3055E-06 
0.2448E-06  0.1997E-06  0.1653E-06  0.1385E-06  0.1174E-06 
0.1004E-06  0.8657E-07  0.7520E-07  0.6574E-07  0.5781E-07 
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0.51 10E-07  0.4537E-07  0.4046E-07  0.2400E-07  0.1514E-07 
0.9963E-08  0.6749E-08  0.4670E-08  0.3283E-08  0.2336E-08 
0.1677E-08  0.1213E-08  0.8824E-09  0.6449E-09  0.4730E-09 
0.3480E-09  0.2568E-09  0.1898E-09 
END 
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